Weak Long-Ranged Casimir Attraction in Colloidal Crystals 
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We investigate the influence of geometric confinement on the free energy of an idealized model 
for charge-stabilized colloidal suspensions. The mean-field Poisson-Boltzmann formulation for this 
system predicts pure repulsion among macroionic colloidal spheres. Fluctuations in the simple 
ions' distribution provide a mechanism for the macroions to attract each other at large separa- 
tions. Although this Casimir interaction is long-ranged, it is too weak to influence colloidal crystals' 
dynamics. 



Experimental evidence collected over 20 years M sug- 
gests that similarly-charged colloidal spheres dispersed in 
water need not simply repel each other. Under some cir- 
cumstances they instead experience an unexpected long- 
ranged attraction. For example, like-charge attractions 
are implicated in the cohesion of metastable superheated 
colloidal crystals [gj3| even though isolated pairs of the 
constituent spheres are observed to repel each other ||,|| . 
Comparable attractions have been measured for pairs 
of spheres confined by two g|| charged planar walls. 
Recent calculations [|7||| reveal that such confinement- 
induced attractions cannot be accounted for by local den- 
sity theory nor by electrohydrodynamic coupling P-pJ ■ 
Such anomalous effects in charge-stabilized colloid there- 
fore challenge our general understanding of interactions 
and dynamics in macroionic systems. 

This Letter addresses fluctuations' contribution to the 
free energy of highly charged colloidal spheres surrounded 
by a neutralizing cloud of small singly-charged countc- 
rions. Highly symmetric monopolar fluctuations in the 
counterion distribution increase the system's free energy. 
We demonstrate that their suppression by boundary con- 
ditions at the spheres' surfaces introduces a long-range 
attraction into the crystal's free energy analogous to the 
Casimir force in quantum electrodynamics, but that it is 
too weak to account for anomalous behavior in charge- 
stabilized suspensions. 

Our treatment is based on the Wigner-Seitz cell model 
which has been studied extensively (TTJ both theoretically 
and through Monte Carlo simulation. It consists of a sin- 
gle spherical macroion of radius a carrying a uniformly 
distributed surface charge — Ze and surrounded by a ther- 
mal cloud of Z point-like counterions at temperature T, 
each carrying a single charge e. The macroion and coun- 
terions are confined by a concentric conducting spherical 
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shell of radius R. This outer shell plays a role analogous 
to the Wigner-Seitz cell boundary in a colloidal crystal. 
More generally, it models the crowding or geometric con- 
finement characteristic of colloidal crystals |Q . 

Previous investigations of this and related models 
|ll| , [l3| have found short-ranged correlation-driven attrac- 
tions between the bounding surfaces under some condi- 
tions, particularly when the counterions are polyvalent. 
They have not found evidence for long-ranged attractions 
in monovalent electrolytes p2[ . 

Our method for evaluating the counterions' partition 
function allows us to investigate much higher macroion 
charges than have been considered before. The outer 
boundary's suppression of counterion fluctuations in- 
duces a long-ranged Casimir-like attraction |l3| between 
the macroion and its neighbors across the Wigner-Seitz 
cell boundary. Although this cell model is far too simple 
to describe the behavior of real charge-stabilized suspen- 
sions, it highlights a previously unexplored mechanism 
for long-ranged confinement-induced like-charge colloidal 
attractions. 

We adopt the path integral formalism reviewed in E3] 
and write the counterions' canonical partition function 
as a functional integral over all possible counterion dis- 
tributions, n(r): 



Q = V2^ZZ Z f er f3f[n] Dn, 



(1) 



where /3 _1 = fcgT is the thermal energy scale at temper- 
ature T, the prime indicates that the number of charges 
is conserved (J n dV = Z) , and Q ] 



f[n] = U[n] + ksT / n \nndV. 



The potential energy functional 
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(2) 



(3) 



describes the counterions' interaction with the local elec- 
tric potential 4>(r). The system's Helmholtz free energy 
is then F= -k B T\nQ. 



One ionic distribution, no(f), minimizes /, and thus 
has the greatest statistical weight in Q. We factor the 
partition function Q = QoQfi into the saddle point con- 
tribution 



Qo 



z z < 
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(4) 



where /q = /[no], and a term Qfi accounting for fluc- 
tuations, Sn, away from uq. Expressing Qfi as a series 
expansion in Sn yields as the lowest-order non-vanishing 
term 



Q f i ~ V2^Z J e- (i5 ' 2f Dn, 



(5) 



where 5 2 f is the second-order change in f[n] due to Sn. 
Terminating this expansion at Gaussian order is justi- 
fied if corrections at higher order in Sn contribute neg- 
ligibly to Qfi. This condition is met if Z is large pa] 
and n (r) itself changes negligibly over the mean radial 
counterion separation: 



1 dn . 2 
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no dr 



(6) 



and can be tested a posteriori once no(r) is evaluated. 
Expanding around any other distribution [|16|,|l7J would 
not give such a criterion for establishing convergence. 

It has been shown [ |l5|Jl7|| that the saddle point corre- 
sponds to the mean-field result 



n = n s e 



-0e<f> 



(7) 



which, combined with the Poisson equation 

V 2 ^--, (8) 

e 

yields the familiar Poisson-Boltzmann (PB) equation, 



-pet 



(9) 



Selecting n s = 3Z/(47ra 3 ) sets the potential's reference 
point conveniently without loss of generality. Eq. (g) 
accounts for the solvent's influence in the so-called prim- 
itive model through its dielectric constant e. 

Following conventional practice, we introduce the Bjer- 
rum length \b = /3e 2 /(47re) and an effective screening 
length k _1 = (47rn s As) -1 / 2 . Taking the system's radial 
symmetry into account leads to 



dr 2 



2d(f> 
r dr 
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-0e<j> 



(10) 



We solve Eq. ( |10| ) subject to two boundary conditions: 
Gauss' theorem at the macroion's surface gives 
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3/?e' 



(11) 
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FIG. 1. Dependence of the mean-field free energy on the 
confining boundary's radius, R, from Eq.(14) with Z = 7300, 
a = 0.326/^m and As = 0.7 nm. The solid curve shows the 
solution of the full PB equation while the dashed line is from 
the linearized PB theory. 



and electroneutrality requires 



dr 



(12) 



Equations ([10|— |l2[) can be solved numerically for cf>(r) 
from which we can calculate the mean-field free energy, 
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r dr 



(14) 

Figure [j] shows /o as a function of R/a, calculated nu- 
merically from the solution of the full PB equation. Since 
fo(R) decreases monotonically with R, the mean-field 
theory predicts pure repulsion. 

An equivalent result obtained by solving the linearized 
PB equation for the same system does have a local min- 
imum for R slightly bigger than a. In the context of col- 
loidal interactions, the linearized mean-field description 
yields the electrostatic component of the pair potential 
due to Derjaguin, Verwey, Landau and Overbeek (DLVO) 
fll8| . The local minimum in the free energy suggests an 
effective electrostatic attraction between like-charged col- 
loidal spheres within the DLVO theory. Comparison with 
the full calculation, however, shows this to be an artifact 
of linearization. 

In addition to confirming the absence of attractions 
in the mean-field description, our numerical results also 
satisfy the condition in Eq. ffl). Thus, we are justified in 
using Eq. (EJ) to calculate Qfi. The second-order change 
in U[n] due to 5n is 



52jJ = ir- I d ^ 
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(15) 



So, from Eq. (§), 
Q/i = yfifcZ J exp 



-/3(5 2 C7 



(fin) 

2n 
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except Ajv.jv which is 1. Separating mode xm from the 
other N — 1 modes in this way allows us to remove the 
explicit constraint: 



We evaluate the functional integral by partitioning the 
system into N concentric shells, 
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where V* is the volume of cell z, rij the number density 
of counterions therein, and n i the equivalent number 
density for the mean-field case. This highly symmet- 
ric partition of the system is appropriate for radial or 
monopolar density fluctuations. A more complete de- 
scription including multipole fluctuations is not neces- 
sary for our purposes because recent calculations reveal 
that these yield only short-ranged attractions (19] . Cross 
terms cancel at the Gaussian level of approximation, so 
we can consider monopole and multipole contributions 
separately. 

We define x, = Srii/ y / 2n 0i /Vi, and divide the system 
so that every cell has the same number of counterions in 
the mean-field distribution, i.e. n^iVi — Z/N . Rescaled 
in this way, 
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Then, numbering the cells from the center outwards so 
that Ti < Tj if % < j, we obtain 
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The prime indicates that the integral is constrained to 
maintain electroneutrality: y'., ccj = 0. Because of this 
condition, we can subtract a constant from each element 
of A without changing the value of A • x. Let us subtract 
Z Xb j (N R) from each element of A to obtain matrix A. 
All of the elements in the last row and column of A vanish 



where x is now a vector of dimension A — 1 , and B is a 
matrix of dimension N — 1 obtained by adding 1 to each 
element in the first A— 1 rows and columns of matrix A. 
B may be expressed as the sum of two matrices, C and 
D, whose components are C^ = 1 + 5, b j and 



_ Z ( 1 1 

B ^~N Xb \--r 



(24) 



where p is the greater of i and j. det C = N, so that 

detB = 7VdetI + CT 1 D, (25) 



where I is the identity matrix. 

Evaluating detB is greatly facilitated if the compo- 
nents of C _1 D are much smaller than 1. C/- 1 = 5ij — l/N 
differs little from the identity matrix. The components 
of D, on the other hand, are bounded above by Dy < 
(Z/N) (X B /a). We previously assumed Z/N > 1 in de- 
riving Eqs. (|J) and (||). But As/a <C 1 for the micron- 
sized spheres in experimental observations, so that we 
may reasonably assume Dy < 1. Even if this were not 
the case, we would be justified in formally taking the 
limit Z/N <§C 1 at this point because the final result can- 
not depend on A. Consequently, 

detB = 7Vexp[Trln(l + Cr 1 D)] (26) 

w Aexp(TrC _1 D). (27) 

In this approximation, the fluctuation contribution to the 
free energy is 

SF= -fc B rTrC _1 D 
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Rewriting the sums over shell indices as integrals over 
radii, we obtain 
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Most of the boundaries, rj, between cells in the mean- 
field distribution are clustered near a. Consequently, 



SF «2vrfc B T / n {r)w (r) 
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r 2 dr, (31) 



where wo(r) = (Att/Z) J no(r')r' dr' is the fraction of 
countcrions within radius r, in the mean field approxi- 
mation. Integrating by parts then yields 
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k B T X B 



l-wg(r) 



dr. 



(32) 



Unlike fo(R), 5F(R) decreases with decreasing R be- 
cause the outer boundary condition suppresses fluctua- 
tions as R approaches a. The resulting attraction there- 
fore is reminiscent of the Casimir attractions previously 
identified in confined electrolytes as well as other systems 
|L3[ . It is interesting to note that monopolar fluctuations 
do not yield an attractive contribution in all geometries; 
the second term in the exponent of Eq. (^8() vanishes 
for unbounded systems such as parallel plates and con- 
centric cylinders. This is consistent with the absence of 
long-ranged like-charge attractions in measurement s P0[ , 
theoretical treatments JlllJlq] , and simulations ]ll| , |l3j ] of 
unbounded systems. Such attractions, therefore, are pe- 
culiar to closed systems, such as the Wigner-Seitz cells of 
colloidal crystals. If there is an R at which the attractive 
force 
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has larger absolute value than the mean-field repulsive 
force 
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(34) 



then F(R) — fo(R) + SF(R) would have a minimum 
at that R. Such a minimum would correspond to a 
fluctuation-mediated bound state for a colloidal crystal of 
nearest neighbor spacing 2R. However, numerical results 
shown in Fig. |] reveal |x /Xr| < 10~ 4 over the entire 
range of conditions studied experimentally. On this basis, 
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FIG. 2. Numerically evaluated ratio \xa/xA maximised 
over R for various values of the macroion radius a and the 
macroion charge Z. A smooth surface is fitted to the data 
points to aid the eye. The plot clearly demonstrates that the 
attraction is four orders of magnitude too small to overcome 
the repulsion. 



we conclude that monopolar fluctuations are not respon- 
sible for the strong and long-ranged attractions reported 
in measurements on charge-stabilized colloid |2|,||,|]|6]] . 

To understand why the ratio |Xa/Xr| is small we esti- 
mate it analytically by assuming that most of the counte- 
rions are clustered close to the macroion surface and that 
only a small fraction p(Z,a,R) of counterions are affected 
by a change in the radius R of the confining shell, and 
that this fraction is uniformly distributed in the volume. 
Evaluating \a (R) to O (p) then yields 



\Xa\ 



Zk B TX B 3p x{x + 2) 



2 2 a 2 (x 2 + l + x) 2 

where x = R/a. Similarly evaluating Xr(R) yields 
Zk B T 3p x 2 



\Xr\ 



2a x 3 -1 



(35) 



(36) 



Truncating to O (p) is justified by numerical investigation 
of the mean-field solution which indicates p < 0.05 in the 
region of interest. In this approximation the ratio 



\Xa\ = X B (x 2 +x-2) l 
\Xr\ a (x 2 + x + 1) x 



(37) 



is independent of p and has a maximium value 0.286 X B /a 
at R = 2.067a. The location and magnitude of the max- 
imum value agrees with the full numerical solution to 
within factors of 1.2 and 4, respectively. This value of 
R would correspond to half the nearest-neighbor separa- 
tion in a fluctuation-stabilized colloidal crystal; it also is 
consistent with the lattice constants observed in super- 
heated metastable colloidal crystallites 0,0] ■ However, 
water has a dielectric constant e = 80 which corresponds 
to a Bjerrum length X B = 0.7 nm at room temperature. 



Consequently, \xa/xA ^ 1 f° r au reasonable colloidal 
radii. Changing solvents would not affect this conclusion 
because of the limited range of accessible values of e. 

We also investigated the possible influence of multiva- 
lent counterions carrying charge qe. Ignoring the rela- 
tively weak Z dependence of Xa/Xr yields Xa/Xr ~ q 2 ■ 
Despite the relative strength of the attraction being 
larger for higher valency counterions, the effect is still 
too weak to induce measurable attractions for physically 
plausible values of q. 

We have demonstrated that suppression of monopo- 
lar ionic fluctuations by surfaces induces a long ranged 
attraction remniscent of Casimir attractions |13|| . This 
interaction is distinct from and complementary to at- 
tractions driven by multipolar fluctuations which have 
been studied elsewhere fL9f . Neither mechanism, how- 
ever, accounts for the strong and long ranged attractions 
observed experimentally between highly charged colloidal 
spheres. 

The long-ranged like-charge attractions observed in 
confined colloid arc not consistent with mean field theo- 
ries for electrolyte structure. Possible explanations must 
incorporate mechanisms such as fluctuations and high- 
order correlations not captured by mean field theory. 
While multipole fluctuations in the distribution of simple 
ions induce strong attractions [ |19| , they are short ranged. 
The present study demonstrates that radially symmetric 
fluctuations can induce long-ranged attractions, but that 
they are too weak to influence colloidal behavior. Conse- 
quently, the explanation must lie in another mechanism 
not yet considered and thus remains an important out- 
standing challenge. 
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